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ABSTRACT 


To  investigate  the  role  of  planetary  beta  on  classical  as  well  as  unique  features  in 
the  northern  Canary  Current  System  (NCCS),  four  numerical  experiments  are  conducted 
with  varying  Coriolis  parameterizations  (/-plane  or  (3-plane).  The  first  two  experiments 
use  a  closed  boundary  and  annual  salinity  forcing  for  the  Mediterranean  Outlfow  (MO). 
The  latter  two  experiments  use  an  open  Mediterranean  Sea  at  the  Strait  of  Gibraltar  and 
seasonal  forcing  for  MO  to  permit  a  more  accurate  investigation  of  the  role  of  the  beta 
effect  on  subsurface  spreading  of  MO  and  Meddies.  All  four  experiments  use  seasonal 
climatological  winds  and  seasonal  thermohaline  gradients  along  the  western  boundary  to 
force  the  model.  Experiments  run  on  a  (3 -plane  (Experiments  2  and  4)  accurately  portray 
classical  eastern  boundary  current  (EBC)  mesoscale  features.  In  addition,  these 
experiments  depict  unique  NCCS  features  associated  with  a  large  embayment  (the  Gulf  of 
Cadiz),  poleward  spreading  of  MO,  and  the  generation  of  Meddies.  Experiments  run  on 
an  /-plane  (Experiments  1  and  3)  show  the  unrealistic  dominance  of  a  continuously 
strengthening  equatorward  jet  that  inhibits  development  of  classical  EBC  and  unique 
NCCS  features.  The  complex  upper  layer  and  subsurface  flow  regimes  of  Experiment  4 
most  realistically  portray  currents,  mesoscale  features  and  Meddies  similar  to  NCCS 
observations. 
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I.  INTRODUCTION 


The  Canary  Current  System  (CCS)  on  the  eastern  boundary  of  the  central  North 
Atlantic  is  a  classical  eastern  boundary  current  (EBC)  system.  Stretching  from  ~10°N  to 
~45°N  along  the  coasts  of  northwest  Africa  and  the  Iberian  Peninsula  (IP),  it  marks  the 
closing  eastern  boundary  of  the  North  Atlantic  Gyre.  Typical  of  other  EBCs,  the  mean 
equatorward  Canary  Current  (CC),  is  a  broad  (-1000  km),  relatively  slow  (-10-30  cm/s), 
yearlong  surface  flow  extending  to  depths  of  -500  m  (Wooster  et  al.,  1976).  The  portion 
of  the  CC  that  stretches  along  the  coast  of  the  IP  is  often  referred  to  as  the  Portugal 
Current  (e.g.,  Tomczak  and  Godfrey,  1994). 

As  in  other  EBCs,  a  poleward  undercurrent  exists  near  the  coast  beneath  the  CC 
(e.g.,  Meincke  et  al.,  1975;  Fiuza,  1980)  as  a  relatively  narrow  (-10-40  km)  and  weak 
(-2-10  cm/s)  flow,  which  is  strongest  between  -100  and  -600  m  depth.  The  depth  and 
strength  of  the  undercurrent  varies  seasonally  and  latitudinally.  In  winter  it  shoals  to  the 
north  near  Cabo  da  Roca  (see  Figures  la  and  lb  for  geographic  locations,  and  Figure  22.1 
of  Barton,  1998  for  bathymetric  contours  and  coastline  geometry  for  the  region),  and 
forms  a  third  flow  component  commonly  referred  to  as  the  Iberian  Current  (IC)  (Haynes 
and  Barton,  1990).  A  narrow  (25-40  km),  relatively  weak  (-20-30  cm/s),  seasonal 
surface  current,  the  IC  is  found  trapped  near  the  coast  against  the  shelf  break  (Fiuza, 
1980;  Frouin  et  al.,  1990;  Haynes  and  Barton,  1990)  and  can  occasionally  be  seen  as  far 
south  as  Cabo  de  Sao  Vincente  (Batteen  et  al.,  2000). 
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A  unique  feature  that  distinguishes  the  CCS  from  other  EBCs  is  the  existence  of 
Mediterranean  Outflow  (MO)  through  the  Strait  of  Gibraltar  into  the  adjacent  Gulf  of 
Cadiz.  A  large  embayment,  the  Gulf  of  Cadiz's  pronounced  east-west  coastline 
orientation  results  in  weaker  upwelling  in  the  Gulf  of  Cadiz  than  to  the  north  or  south  of 
the  Gulf  of  Cadiz,  due  to  the  dominant  equatorward  trade  wind  direction.  The  Gulf  of 
Cadiz  also  creates  a  large  separation  between  the  two  west  coast  upwelling  regimes  so 
that  no  continuous  flow  between  the  two  appears  to  exist  (Barton,  1998). 

The  salty  MO  plume  that  exits  the  Mediterranean  Sea  through  the  Strait  of 
Gibraltar  is  diluted,  thickens,  and  becomes  vertically  differentiated  into  two  distinct  cores 
as  it  flows  westward  into  the  Gulf  of  Cadiz  (Iorga  and  Lozier,  1999).  At  ~7°W  in  the 
Gulf  of  Cadiz,  both  a  shallow  core  at  depths  of  -600-900  m  and  a  deeper  core  at  —1100- 
1200  m  exist.  Both  cores  continue  to  flow  westward  along  the  southern  coast  of  Spain 
and  turn  poleward  around  Cabo  de  Sao  Vincente  (Amber  and  Howe,  1979;  Iorga  and 
Lozier,  1999).  A  third,  shallower,  poleward  core  of  Mediterranean  water  has  also  been 
traced  from  the  Strait  of  Gibraltar  northward  to  ~38.5°N  off  western  Portugal  (Amber, 
1982).  In  addition,  a  climatological  cyclonic  circulation  in  the  southwestern  Gulf  of 
Cadiz  acts  to  spread  salty  MO  south  of  ~34°N  (Iorga  and  Lozier,  1999). 

The  northern  CCS  (NCCS)  is  influenced  predominantly  by  equatorward, 
upwelling  favorable  winds  produced  by  the  eastern  half  of  the  Azores  High.  Located  in 
the  northeastern  part  of  the  Atlantic  Ocean,  the  Azores  High  is  a  semi-permanent 
subtropical  high  pressure  system  similar  in  nature  and  behavior  to  the  North  Pacific 
Subtropical  High  (Nelson,  1977).  As  such,  the  center  of  the  Azores  High  migrates 
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meridionally  with  the  seasons,  reaching  its  southernmost  extent  near  27°N  in  March  and 
ridging  north  to  ~33°N  by  August.  The  building  and  migration  of  the  Azores  High  cause 
wind  stress  values  over  the  NCCS  region  to  vary  temporally,  resulting  in  both  alongshore 
and  cross-shore  variability  of  the  predominantly  equatorward  winds  along  the  west  coast 
of  the  IP  and  the  northwest  coast  of  Africa.  The  east-west  pressure  contrast  between 
Portugal  and  the  center  of  the  Azores  High  during  summer  is  ~8  mb.  In  winter,  this 
pressure  gradient  weakens  to  ~1  mb.  As  a  result  of  this  seasonally  changing  pressure 
gradient,  considerably  stronger  northerly  and  northwesterly  winds  occur  in  the  NCCS 
region  in  summer,  while  northerly  and  northwesterly  winds  become  weaker  to  even  slight 
southerly  off  the  Iberian  Peninsula  in  winter  (Batteen  et  al.,  2000).  The  shift  in 
maximum  wind  stress  also  causes  upwelling  favorable  winds  to  shift  from  ~27°N  near 
the  Canary  Islands  in  January,  to  ~43°N  off  Portugal  by  July  (Fiuza,  1982). 

Like  other  classical  EBCs,  observations  of  the  sea  surface  in  the  NCCS  region 
have  shown  highly  energetic  mesoscale  features  such  as  jet-like  surface  currents, 
meanders,  eddies  and  filaments  over  the  broad  climatological  mean  flow  of  the  CC. 
Satellite  sea  surface  images  have  shown  nearshore  upwelling  during  periods  of  upwelling 
favorable  winds  with  several  narrow  filaments  of  cooler  water  extending  off  the  coast  of 
the  Iberian  Peninsula  (Fiuza  and  Sousa,  1989)  and  Cape  Ghir  in  northwest  Africa  (Van 
Camp  et  al.,  1991;  Hagen  et  al.,  1996).  In  these  images,  upwelling  filaments  often 
extend  -80-150  km  offshore,  with  alongshore  spacing  of  -80-100  km  between  filaments 
(Haynes  et  al.,  1993),  and  often  terminate  with  dipole  eddy  pairs  (Fiuza  et  al.,  1982; 
Barton,  1998).  Observations  have  also  shown  anticyclonic  and  cyclonic  pairs  of 
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mesoscale  eddies  on  the  order  of  100  km  off  the  IP  coast  (Fiuza,  1984;  Stammer  et  al., 
1991).  These  mesoscale  features  have  been  observed  during  periods  of  predominantly 
upwelling  favorable  winds  and  appear  to  be  located  near  prominent  coastline 
irregularities  such  as  capes.  These  observations  provide  evidence  that  wind  forcing  along 
coastline  irregularities  appear  to  be  important  mechanisms  in  the  formation  and 
sustainment  of  many  of  the  mesoscale  features  found  in  the  NCCS  domain  and  other 
EBC  regions  (Batteen  et  al,  2000). 

Unique  to  the  NCCS  is  the  generation  of  anticyclonic  submesoscale  coherent 
vortices  (SCVs)  or  Meddies.  Numerical  studies  suggest  that  baroclinic  instability  of  the 
northward  dense  plume  of  salty  MO  along  the  IP  continental  slope  leads  to  the  generation 
of  Meddies  (Kase  et  al.,  1989).  As  a  result  of  numerous  observations  over  the  past 
decade,  the  primary  generation  region  of  Meddies  is  widely  accepted  to  be  near  Cabo  de 
Sao  Vincente,  off  southwest  Portugal.  Several  different  trajectories  of  Meddies  have 
been  observed,  including  a  southwestward  movement  into  the  Canary  Basin,  and 
westward  translations  south  of  the  Azores  (Richardson  and  Tychensky,  1998). 

Over  the  last  two  decades  numerous  modeling  studies  have  focused  on  the  driving 
mechanisms  of  complex  mesoscale  activities  in  EBC  regions,  including  upwelling 
filaments,  highly  energetic  eddies,  and  meandering  jets.  For  the  California  Current 
System,  the  classical  EBC  system,  Ikeda  et  al.  (1984a,  b)  and  Haidvogel  et  al.  (1991) 
studied  baroclinic  instability,  coastline  irregularities,  and  bottom  topography  as  possible 
mechanisms,  while  Batteen  et  al.  (1989),  McCreary  et  al.  (1991),  Pares-Sierra  et  al. 
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(1993),  Batteen  (1997)  and  Batteen  et  al.  (2000)  studied  wind  forcing  as  a  possible 
generation  mechanism. 

For  over  a  decade,  a  high  resolution,  multi-level,  primitive  equation  model  has 
been  used  to  examine  the  response  to  wind  forcing  of  an  idealized  flat-bottomed  oceanic 
regime  along  different  eastern  ocean  boundaries  (e.g.,  Batteen,  1989;  Batteen,  1997, 
Batteen  et  al.,  2000).  Results  from  sensitivity  studies  have  produced  various  results 
depending  on  the  type  of  wind  forcing,  coastline  geometry  and  Coriolis  parameterization 
used.  A  summary  of  experimental  design  for  each  of  these  studies  is  contained  in  Tables 
1, 2  and  3. 

For  the  California  Current  System,  Batteen  (1989)  used  a  band  of  steady 
alongshore,  upwelling  favorable  winds,  either  with  or  without  alongshore  variability,  as 
forcing  on  both  an /-plane  and  (3-plane  (see  Table  1).  This  study  was  conducted  using  a 
straight  California  coast  only.  Eddies  and  jets  with  strong  onshore  and  offshore  directed 
flows  were  generated,  along  with  a  coastal  jet  and  undercurrent  on  a  scale  and  magnitude 
consistent  with  available  observed  data.  The  results  supported  the  hypothesis  that 
alongshore  varying  wind  forcing  on  a  (3 -plane  is  an  important  generation  mechanism  for 
eddies  in  the  California  Current  System. 

For  the  California  Current  System,  Batteen  (1997)  forced  the  model  from  rest 
using  seasonal  climatological  winds  along  a  straight  California  coast  on  both  an  /-plane 
and  (3-plane  (see  Table  2).  This  study  also  included  seasonal  climatological  wind  forcing 
on  an  irregular  coastline  using  (3-plane  only  and  focused  on  the  role  of  the  beta  effect  in 
generating,  evolving  and  maintaining  classical  mesoscale  features  in  the  California 
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Current  System.  The  results  of  this  study  showed  that  a  surface  equatorward  current, 
upwelling  of  cooler  water  along  the  coast,  and  a  poleward  undercurrent  develop  in 
response  to  the  prevailing  wind  direction.  In  addition,  baroclinic  and  barotropic 
instabilities  develop  in  the  equatorward  surface  current  and  poleward  flow,  resulting  in 
the  generation  of  meanders  near  the  coast,  upwelling  filaments  extending  offshore,  and 
the  subsequent  development  of  seasonal  or  semi-permanent  cyclonic  and  anticyclonic 
eddies  farther  offshore.  While  the  meridional  variability  of/(3-plane)  and  variability  in 
the  equatorward  component  of  wind  stress  were  shown  to  be  necessary  ingredients  for 
generating  realistic  vertical  and  horizontal  structures  for  the  cores  of  the  surface 
equatorward  and  subsurface  poleward  currents,  irregularities  in  the  coastline  geometry 
were  shown  to  be  important  for  “anchoring”  upwelling  and  filaments,  and  for  enhancing 
the  growth  of  meanders  and  eddies.  The  results  from  several  numerical  experiments  in 
this  process-oriented  study  supported  the  hypothesis  that  wind  forcing  and  coastline 
irregularities  on  a  3-plane  are  important  mechanisms  for  generation  of  observed 
mesoscale  features  of  the  California  Current  System  (Batteen,  1997). 

In  a  more  recent  study  Batteen  et  al.,  (2000)  used  the  results  of  four  numerical 
experiments  of  increasing  complexity  to  investigate  the  unique  and  classical  features  of 
the  NCCS  (see  Table  3).  All  four  experiments  were  run  on  a  3_pl&ne.  Experiment  1 
included  seasonal  wind  forcing  on  a  straight  coast,  while  Experiments  2  and  3  added  the 
effects  of  irregular  coastline  geometry  to  the  seasonal  wind  forcing.  The  most  complex 
experiment.  Experiment  4,  added  the  effects  of  both  thermohaline  gradients  and  annual 
MO  to  the  seasonal  wind  forcing  and  irregular  coastline  geometry.  Consistent  with 
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previous  studies,  wind  forcing  was  shown  to  be  a  key  generative  mechanism  for  EBC 
mesoscale  features,  while  capes  were  shown  to  be  key  areas  for  enhanced  upwelling, 
extensive  filaments,  maximum  current  velocities  and  enhanced  growth  of  cyclonic 
meanders  and  eddies.  Unique  features  of  the  NCCS,  including  the  development  of 
anticyclonic  meanders  and  eddies  in  an  embayment  like  the  Gulf  of  Cadiz  and  the 
generation  of  Meddies  from  MO,  were  also  highlighted  (Batteen  et  al,  2000). 

The  objective  of  this  study  is  to  build  on  previous  studies  and  to  isolate  the 
dependence  of  Coriolis  parameterization  (3-plane  vs.  /-plane)  on  classical  as  well  as 
unique  mesoscale  features  in  the  NCCS.  Specifically,  the  generation,  evolution,  and 
maintenance  of  currents,  meanders,  eddies,  filaments,  MO,  and  Meddies  is  studied.  The 
results  of  several  numerical  experiments  (see  Table  4)  using  a  high-resolution,  multi-level 
PE  model  are  explored.  In  all  experiments,  irregular  coastline  geometry  is  used  and  the 
model  is  forced  using  seasonal  climatological  winds  with  both  alongshore  and  cross¬ 
shore  variability.  In  addition,  each  experiment  includes  the  effects  of  thermohaline 
forcing  on  the  western  boundary  and  either  annual  or  seasonal  salinity  forcing  of  MO. 
Experiments  1  and  2  are  forced  using  annual  values  of  salinity  near  the  Strait  of  Gibraltar 
to  include  the  effects  of  MO  on  the  NCCS.  In  both  experiments,  the  eastern  Gulf  of 
Cadiz  exists  as  a  closed  boundary.  Experiment  1  introduces  the  effect  of  running  the 
model  on  an /-plane  for  the  NCCS.  Experiment  2  is  the  same  as  Experiment  1  except  the 
model  is  run  on  a  (3-plane.  The  results  of  Experiment  2  are  compared  to  the  results  of 
Experiment  1  to  explore  the  role  of  the  planetary  beta  effect  on  the  NCCS.  In  order  to 
more  realistically  examine  the  effects  of  MO  on  the  NCCS,  Experiments  3  and  4  are 
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forced  using  seasonal  values  of  salinity  for  the  MO.  In  addition,  the  Strait  of  Gibraltar  is 
opened  up  to  allow  for  the  free  exchange  of  water  in  and  out  of  the  Mediterranean  Sea. 
Experiment  3  investigates  the  effect  of  MO  when  the  model  is  run  on  an /-plane,  while 
Experiment  4  is  the  same  as  Experiment  3  except  the  model  is  run  on  a  3-plane.  The 
results  of  Experiment  4  are  compared  to  the  results  of  Experiment  3  to  explore  the  role  of 
planetary  beta  effect  on  MO  and  Meddies  in  the  NCCS.  The  results  of  Experiment  4,  the 
most  realistic  experiment,  are  also  compared  to  observations  of  surface  and  sub-surface 
mesoscale  activities  in  the  NCCS. 

This  study  is  organized  as  follows.  In  Chapter  II  we  describe  the  numerical 
model  and  the  specific  experimental  conditions.  The  results  of  the  numerical 
experiments  are  presented  in  Chapter  HI  along  with  a  comparison  with  observations.  A 
summary  is  presented  in  Chapter  IV. 
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II.  MODEL  DESCRIPTION 


A.  MODEL  EQUATIONS 

This  study  uses  the  eddy-resolving,  limited-area  EBC  model  of  Batteen  (1997) 
with  open  northern,  western,  and  southern  boundaries.  Additionally,  for  this  study,  the 
model  has  been  adapted  to  have  an  open  eastern  boundary  for  the  Strait  of  Gibraltar. 

The  model  is  multi-level,  uses  non-adiabatic  primitive  equations  on  a  (3-plane  or 
/-plane,  and  has  both  baroclinic  and  barotropic  velocity  components.  It  is  based  on  the 
hydrostatic,  Boussinesq,  and  rigid  lid  approximations.  Equations  governing  the  model 
are  as  follows: 


du  -Idp  ,  „4  „  d2u 

dt~p0dx+fV  A**u  +  K"&* 

(1) 

dv  -Idp  4  d2v 

dt~p0dy  **  4  v  + 

(2) 

du  dv  dw  ^ 

!k+!fy+!k =0 

(3) 

dp 

(4) 

p  =  p0ll-a(T-T„)+^S-S<in 

(5) 

S.  _A  v't+js:  }2T 

dt~  A«yT  +  K^dz2 

(6) 

dS  4  d2S 

—  =  -Ay<S  +  KH^T 
dt  H  H  dz2 

(7) 
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In  the  above  equations,  t  is  time,  T  is  temperature,  S  is  salinity,  p  is  density,  and  p  is 
pressure.  A  right-handed  Cartesian  coordinate  system  (x,y,z)  is  used  where  x  points 
toward  shore,  y  alongshore,  and  z  upward.  The  corresponding  velocity  components  are 
(«,v,w).  Table  5  provides  a  list  of  symbols  found  in  the  model  equations,  as  well  as 
values  of  constants  used  throughout  the  study. 

A  space-staggered  B-scheme  (Arakawa  and  Lamb,  1977)  is  used  for  the 
horizontal  finite  differencing.  Batteen  and  Han  (1981)  have  shown  this  scheme  to  be 
appropriate  when  the  grid  spacing  is  approximately  on  the  same  order  as,  or  less  than,  the 
Rossby  radius  of  deformation.  The  horizontal  grid  spacing  is  9  km  in  the  east-west 
direction  and  11  km  in  the  north-south  direction,  while  the  internal  Rossby  radius  of 
deformation  is  ~30  km,  as  determined  by  Feliks  (1985).  The  model  uses  ten  vertical 
layers,  with  constant  z-levels,  at  depths  of  10, 30, 75, 150, 250, 400,  600,  1226,  2283,  and 
3656  m.  This  vertical  spacing  scheme  of  the  model  concentrates  more  layers  above  the 
thermocline,  which  is  in  the  dynamically  active  portion  of  the  ocean. 

The  model  domain  (Figure  la)  encompasses  the  west  coasts  of  Iberia  and 
Morocco,  from  30°  N  to  42.5°  N  (1408  km  alongshore),  and  from  5°  W  to  17.5°  W  (1 152 
km  cross-shore).  The  coastal  boundaries  of  the  model  domain  are  closed  for  the  first  two 
experiments  (e.g.,  see  Figure  la)  and  then  opened  at  the  Mediterranean  for  the  second 
two  experiments  (e.g.,  see  Figure  lb).  Where  closed,  the  coastal  boundaries  have  both 
the  tangential  and  normal  components  of  velocity  set  to  zero.  Bottom  topography  has 
been  omitted  to  focus  on  the  roles  of  climatological  wind  forcing,  thermohaline  gradients, 
Coriolis  parameterization,  and  the  MO.  The  eastern  boundary  is  modeled  a  vertical  wall 
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(Figure  la)  or  as  an  open  boundary  at  the  Strait  of  Gibraltar  (Figure  lb)  to  account  for 
more  realistic  exchange  with  the  Mediterranean  (Figure  lb).  The  constant  depth  used  in 
the  model  is  4500  m. 

The  northern,  southern,  and  western  borders,  as  well  as  a  portion  of  the  eastern 
border  in  two  experiments,  are  open  boundaries,  which  use  a  modified  version  of  the 
radiation  boundary  conditions  of  Camerlengo  and  O'Brien  (1980).  Biharmonic  lateral 
heat  and  momentum  diffusion  is  used  in  the  model  with  the  same  choice  of  coefficients 
17  4 

(i.e.,  2.0  xlO  cm  /s)  as  in  Batteen  (1997).  Holland  (1978)  showed  that  highly  scale- 
selective  biharmonic  diffusion  acts  predominantly  on  submesoscales,  while  Holland  and 
Batteen  (1986)  found  that  baroclinic  mesoscale  processes  could  be  damped  by  Laplacian 
lateral  heat  diffusion.  As  a  result,  the  use  of  biharmonic  lateral  diffusion  should  allow 
mesoscale  eddy  generation  via  barotropic  (horizontal  shear)  and/or  baroclinic  (vertical 

shear)  instability  mechanisms.  Weak  (0.5  cm2/s)  vertical  eddy  viscosities  and 
conductivities  are  also  used.  Bottom  stress  is  parameterized  by  a  simplified  quadratic 
drag  law  (Weatherly,  1972),  as  in  Batteen  (1997). 

The  method  of  solution  is  straightforward  with  the  rigid  lid  and  flat  bottom 
assumptions  because  the  vertically  integrated  horizontal  velocity  is  subsequently 
nondivergent.  The  vertical  mean  flow  can  be  described  by  a  streamfunction  which  can  be 
predicted  from  the  vorticity  equation,  while  the  vertical  shear  currents  can  be  predicted 
after  the  vertical  mean  flow  is  subtracted  from  the  original  equations.  The  other 
variables,  i.e.  temperature,  density,  salinity,  vertical  velocity,  and  pressure,  can  be 
explicitly  obtained  from  the  thermodynamic  energy  equation  (6),  equation  of  state  (5), 
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salinity  equation  (7),  continuity  equation  (3),  and  hydrostatic  equation  (4),  respectively 
(For  more  complete  details  on  the  method  of  solution,  see  Batteen,  1997). 

B.  FORCING  CONDITIONS  AND  EXPERIMENTAL  DESIGN 

The  model  is  forced  from  rest  with  climatological  wind  fields  from  a  2.5°  by  2.5° 
grid  of  the  European  Centre  for  Medium  Range  Weather  Forecasts  (ECMWF)  near¬ 
surface  wind  analyses  (Trenberth  et  al.,  1990).  The  monthly  mean  stresses  based  on 
twice  daily  wind  analyses  from  1980-1989  have  been  interpolated  spatially  to  the  9  by  1 1 
km  model  resolution  and  temporally  to  daily  wind  values.  Sample  wind  fields  used  can 
be  seen  in  Figure  2,  which  shows  the  annual  migration  of  the  Azores  Subtropical  High 
from  the  south  in  the  winter  (e.g.,  Figure  2a)  to  its  most  northern  extent  in  the  summer 
(e.g..  Figure  2d).  The  atmospheric  pressure  pattern  for  November  (not  shown)  and 
December  (Figure  2a)  depicts  a  low  to  the  north  and  the  Azores  High  to  the  south.  This 
results  in  a  wind  divergence  with  weakly  poleward  winds  north  of  ~40°N  and 
equatorward  winds  to  the  south.  During  January  (not  shown)  and  February  (Figure  2b) 
the  divergence  in  the  wind  field  migrates  poleward.  By  March  (not  shown)  and  April 
(Figure  2c)  an  equatorward  component  in  the  wind  field  is  observed  along  the  entire 
domain.  The  strongest  equatorward  winds  are  discernible  from  July  (not  shown)  through 
August  (Figure  2d).  In  September  (not  shown)  and  October  (Figure  2e)  the  winds  start  to 
weaken  throughout  the  domain,  and  divergence  in  the  wind  field  is  observed  in  the  north 
by  October  (Figure  2e).  By  November  (not  shown)  the  wind  divergence  returns  to  ~40°N. 
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Table  4  summarizes  the  design  of  the  four  experiments.  Note  that  in  all  the 
experiments  the  model  is  forced  from  rest  with  seasonal  ECMWF  winds.  The  model  is 
then  allowed  to  reach  a  quasi-equilibrium  state  by  running  each  experiment  for  three 
years. 

In  addition,  in  all  the  experiments,  seasonal  temperature  and  salinity 
climatological  conditions  for  the  upper  seven  levels  from  Levitus  and  Boyer  (1994)  and 
Levitus  et  al.  (1994)  are  used  to  initialize  the  model,  and  every  third  day,  to  force  the 
model  at  the  western  boundary  at  17.5°W.  The  seasonal  temperature  and  salinity-forcing 
conditions,  which  are  initially  assumed  to  be  homogeneous,  are  shown  for  one  of  the 
seven  levels  in  Figure  3  for  the  northern  (42.5°N)  and  southern  (30°N)  boundaries  of  the 
model  domain.  Since  the  lower  three  levels  do  not  exhibit  much  horizontal  variation, 
they  are  assumed  to  be  constant  for  each  level.  The  temperature  values  used  for  levels  8- 
10  are  8.52°C,  3.59°C,  and  2.09°C.  The  salinity  used  for  the  lower  three  levels  is  34.7. 

While  the  temperatures  to  the  south  are  warmer  than  those  to  the  north  above 
1000  m  depth,  only  the  upper  level  temperature  conditions  (e.g..  Figure  3)  show 
significant  seasonal  variability  with  a  temperature  maximum  in  September  and  a 
temperature  minimum  in  February  throughout  the  whole  region.  Below  ~150m  depth 
(not  shown),  both  the  seasonal  temperature  fluctuations  and  the  temperature  gradient 
weaken,  as  expected.  In  contrast,  the  salinity  conditions  above  1000  m  (e.g..  Figure  3), 
which  in  the  upper  levels  show  less  (more)  saline  water  to  the  north  (south),  have  no 
significant  seasonal  cycle. 
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In  Experiments  1  and  2,  the  effects  of  the  MO  on  the  NCCS  are  included  while 
keeping  the  eastern  coastal  boundary  closed  (see  Figure  la).  The  salinity  at  each  level  for 
the  MO  is  input  as  a  constant  (annual)  value  for  five  grid  points  at  -36°  N  every  time  step 
with  the  values  given  in  Table  6  (from  Levitus  et  al,  1994). 

As  Table  4  shows,  Experiment  1  examines  the  response  of  seasonal  wind  forcing 
and  the  effects  of  thermohaline  gradients  on  an  irregular  coastline  when  the  model  is  run 
on  an /-plane.  To  isolate  the  effects  of  planetary  beta,  Experiment  2  is  run  using  the  same 
initialization  and  forcing  parameters  as  in  Experiment  1,  but  on  a  (5-plane. 

To  more  realistically  encompass  the  effects  of  MO  on  the  NCCS,  the  eastern 
coastal  boundary  is  opened  at  five  grid  points  to  mask  the  Strait  of  Gibraltar  in 
Experiments  3  and  4  (see  Figure  lb).  Seasonal  values  of  salinity  are  forced  from  the 
upper  five  levels  for  the  MO  every  third  day,  while  salinity  values  for  the  lower  levels  are 
input  as  a  constant.  The  seasonal  values  forced  from  the  upper  five  levels  for  one  year  are 
shown  in  Figure  4  (from  Levitus  et  al,  1994). 

Experiments  3  and  4  are  run  to  isolate  the  effects  of  planetary  beta  on  MO  in  the 
NCCS.  Experiment  3  investigates  the  behavior  of  MO  when  the  model  is  run  on  an 
/-plane,  while  Experiment  4  is  run  on  a  3 -plane. 
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III.  RESULTS  FROM  MODEL  SIMULATIONS 


A.  EXPERIMENT  1  -  EFFECTS  OF /-PLANE  ON  IRREGULAR  COASTLINE 
GEOMETRY 

In  Batteen  (1989)  and  Batteen  (1997),  several  experiments  were  run  on  an /-plane. 
In  all  experiments  straight  coastline  geometry  was  used,  while  the  model  was  forced 
using  different  types  of  wind  stress.  This  experiment,  Experiment  1,  is  also  run  on  an  /- 
plane  but  includes  a  realistic  irregular  coastline  as  well  as  being  forced  by  seasonal 
climatological  winds,  thermohaline  forcing  at  the  western  boundary,  and  annual  salinities 
for  MO. 


1.  Spin-up  Phase 

Initially,  the  combination  of  wind  forcing  and  thermohaline  gradients  results  in 
different  oceanic  responses  depending  on  the  season.  As  the  temperature  decreases 
poleward  (Figure  3),  there  is  an  increase  of  density  poleward,  and  the  corresponding 
pressure  gradient  establishes  an  onshore  geostrophic  inflow  from  the  ocean  interior.  In 
the  winter  (e.g.  Figure  5a),  in  the  poleward  end  of  the  model  domain,  the  large  pressure 
gradient  and  poleward  wind  stress  result  in  onshore  geostrophic  and  Ekman  flow,  and  a 
coastal  poleward  surface  current.  Off  Morocco,  in  the  equatorward  end  of  the  model 
domain,  the  smaller  pressure  gradient  and  equatorward  wind  stress  result  in  weak  onshore 
geostrophic  flow,  offshore  Ekman  flow,  and  a  coastal  equatorward  surface  current. 
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The  combined  effects  of  an  increase  in  equatorward  winds  along  the  entire  coastal 
domain  during  the  peak  upwelling  season  (July  through  September)  and  a  weakened 
pressure  gradient,  lead  to  substantial  strengthening  of  equatorward  flow.  This  flow 
remains  close  to  the  west  coast  of  Iberia  and  Morocco  and  meanders  in  a  slight 
anticyclonic  manner  as  it  bisects  the  Gulf  of  Cadiz  (e.g.  Figure  5b).  A  large  quasi¬ 
permanent  cyclonic  eddy,  extending  to  -1000  m  depth  (e.g.  Figure  5c),  forms  in  the 
eastern  half  of  the  Gulf  of  Cadiz  to  the  north  and  east  of  the  anticyclonic  equatorward 
flow  (e.g.  Figure  5b).  The  weakening  (strengthening)  of  onshore  (offshore)  flow  results 
in  upwelled  water  replacing  the  coastal  waters  along  the  west  coasts  of  Iberia  and  Africa. 
The  coldest  water  is  found  near  coastal  promontories  (e.g.,  off  Cape  Ghir  and  Cape 
Beddouzza  in  Figure  5d). 

A  typical  cross-section  (Figure  6a)  of  the  currents  shows  a  relatively  strong 
coastal  equatorward  jet  extending  to  -1000  m  depth.  The  surface  current  lies  within  -30 
km  of  the  coast,  and  has  core  velocities  of  -35-40  cm/s.  There  is  a  deep  (-2800  m  depth) 
poleward  undercurrent,  also  within  -30  km  of  the  coast,  below  the  equatorward  jet,  with 
a  core  velocity  of  -6  cm/s. 

By  fall,  a  typical  cross-section  of  the  currents  in  the  poleward  end  of  the  model 
domain  north  of  ~41°N  (e.g.  Figure  6b)  shows  that  the  equatorward  coastal  jet  (-30  km 
from  the  coast)  continues  to  maintain  intensity,  with  core  velocities  of  -40  cm/s.  In 
addition,  at  -2200  m  depth,  the  poleward  undercurrent  lies  within  -30  km  of  the  coast 
and  has  a  core  velocity  of  -10  cm/s.  Instead  of  the  expected  return  of  a  surface  poleward 
flow  in  the  northern  portion  of  the  model  domain  in  response  to  the  return  to  poleward 
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wind  stress  in  the  region,  the  surface  equatorward  coastal  jet  continues  to  dominate  the 
entire  model  domain  from  north  to  south. 

A  typical  cross-section  of  salinity  near  Cabo  de  Sao  Vincente  (e.g.  Figure  7) 
shows  a  salinity  maximum  -35.71  at  -1200  m  depth  associated  with  annual  forcing  of 
salinity  for  MO.  A  salinity  minimum  of  35.53,  discernible  at  -600  m  depth,  is  associated 
with  the  equatorward  coastal  jet  (not  shown). 

2.  Quasi-equilibrium  Phase 

A  continuation  of  the  model  run  for  Experiment  1  shows  that  the  system  has 
reached  quasi-steady  state  by  the  end  of  year  2  (not  shown).  Using  years  2  and  3  of  the 
longer  experimental  run,  the  model  shows  further  intensification  of  the  coastal 
equatorward  jet  with  time  and  no  differences  in  oceanic  responses  based  on  the  season. 
The  cyclonic  eddy  in  the  Gulf  of  Cadiz  migrates  and  is  incorporated  into  the  mean 
equatorward  flow,  while  an  anticyclonic  eddy  develops  in  its  wake  (e.g.,  Figure  8a).  As 
the  strength  of  the  equatorward  jet  increases,  the  anticyclonic  meander  of  the 
equatorward  flow  seen  earlier  is  replaced  by  a  more  north-south  orientated  flow  (Figure 
8b).  During  years  2  and  3,  there  are  periodic  events  in  which  the  anticyclonic  eddy  in  the 
Gulf  of  Cadiz  is  forced  south  by  the  equatorward  jet  out  of  the  Gulf  and  into  the  mean 
equatorward  flow  of  the  jet.  During  each  of  these  instances,  a  new  anticyclonic  eddy 
forms  in  the  same  region  of  the  Gulf  of  Cadiz  as  the  anticyclonic  eddy  seen  earlier  (e.g. 
Figure  8c). 
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In  summer,  (e.g.  Figure  9),  upwelling  is  visible  along  the  west  coast  of  Iberia  and 
along  the  coast  of  Morocco,  reaching  a  maximum  in  October.  Areas  of  enhanced 
upwelling  occur  near  prominent  capes  to  the  north  (e.g.,  off  Cabo  da  Roca  in  Figure  9) 
and  off  Morocco. 

A  typical  cross-section  of  the  currents  during  the  fall  and  winter  (e.g.  Figure  10a) 
shows  the  strong  coastal  equatorward  jet  extending  to  -1100  m  depth.  The  surface 
current  is  within  -40  km  of  the  coast  and  has  core  velocities  of  -70-80  cm/s.  The 
poleward  undercurrent  at  -2800  m  depth  is  within  -70  km  of  the  coast,  with  a  core 
velocity  in  excess  of  -10  cm/s.  A  typical  cross  section  of  currents  in  the  poleward  end 
of  the  model  domain  north  of  ~41°N  (e.g..  Figure  10b)  shows  a  surface  equatorward  jet 
close  to  the  shore  instead  of  the  expected  surface  poleward  flow  in  response  to  the 
poleward  wind  stress. 

A  typical  cross-section  of  temperature  during  year  3  (e.g.  Figure  11a),  shows 
temperature  patterns  consistent  with  upwelling  between  -200  and  -1200  m  depths  in 
response  to  the  continuous  and  deep  equatorward  jet.  By  year  3,  the  coastal  jet  acts  as  a 
thermal  boundary,  separating  warm  waters  offshore  of  the  equatorward  jet  from  cooler 
waters  along  the  coast  and  in  the  Gulf  of  Cadiz.  For  example,  Figure  lib  shows  that  a 
typical  temperature  for  October  west  of  the  equatorward  jet  is  -22°  C,  while  a  typical 
temperature  to  the  east  of  the  jet  is  -18°  C  (excluding  coastal  upwelling  areas). 

The  coastal  jet  behaves  as  a  similar  boundary  for  salinity,  separating  fresher 
waters  offshore  of  the  equatorward  jet  from  saltier  waters  along  the  coast  and  in  the  Gulf 
of  Cadiz.  For  example,  Figure  12  shows  that  a  typical  salinity  for  October,  at  600  m 
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depth,  west  of  the  equatorward  jet  is  less  than  -35.55  while  a  typical  salinity  to  the  east 
of  the  jet  is  greater  than  -35.65.  In  addition,  a  typical  cross-section  of  salinity  off  Cabo 
de  Sao  Vincente  (e.g.  Figure  13)  shows  a  salinity  maximum  of  -35.71,  associated  with 
annual  salinity  forcing  for  MO,  inshore  of  the  deep  equatorward  coastal  jet  at  -650  m 
depth,  and  a  salinity  minimum  of  -35.53  offshore  of  the  jet  at  around  the  same  depth. 

3.  Discussion 

Many  of  the  results  of  this  /-plane  experiment  for  the  NCCS  are  consistent  with 
previous /-plane,  seasonally  wind-driven  experiments  (e.g.,  Batteen,  1989  and  Batteen, 
1997),  which  used  straight  coastline  geometry  for  the  California  Current  System.  For 
example,  in  this  experiment  and  in  previous  cases,  the  seasonal  climatological  wind  field 
intensifies  the  surface  equatorward  current  and  deepens  and  weakens  the  poleward 
undercurrent.  As  a  result,  there  are  no  strong  vertical  and  horizontal  shears  in  the  upper 
layers  between  the  surface  equatorward  and  the  subsurface  poleward  currents. 
Consequently,  no  meanders,  filaments  or  eddies  develop  along  the  coast  in  the  vicinity  of 
the  equatorward  jet. 

An  exception  to  the  lack  of  mesoscale  activity  occurs  in  the  Gulf  of  Cadiz,  which 
is  a  large  embayment  and  a  unique  feature  of  the  NCCS.  The  sharp  angle  in  the  coastline 
at  Cabo  de  Sao  Vincente,  the  large  extent  of  the  embayment,  and  the  increasing 
momentum  of  the  deep  equatorward  jet,  cause  the  equatorward  flow  to  bypass  most  of 
the  Gulf  of  Cadiz.  The  Gulf  of  Cadiz  becomes  a  spawning  region  for  predominantly 
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anticyclonic  mesoscale  activity.  (The  cyclonic  eddy  that  forms  in  the  northeastern  part  of 
the  Gulf  of  Cadiz  (e.g.,  Figure  5b)  can  be  attributed  to  the  onshore  pressure  gradient  and 
poleward  flow  that  occurs  before  the  equatorward  jet  intensifies.)  The  series  of 
anticyclonic  eddies  that  periodically  develop  before  being  caught  up  in  the  equatorward 
flow  during  the  quasi-equilibrium  phase  are  due  to  eastward  (westward)  flow  along  the 
northern  (southern)  coasts  of  the  Gulf  and  equatorward  flow  on  the  west  coast  of  the 
Gulf,  resulting  in  formation  of  anticyclonic  eddies  in  the  Gulf  of  Cadiz  (e.g.  Figure  9c). 

The  widths  of  the  coastal  jet  and  undercurrent  are  approximately  confined  to 
within  the  first  internal  Rossby  radius  of  deformation  of  the  coast  along  Iberia  and 
Morocco,  in  concurrence  with  the  expected  behavior  of  eastern  boundary  currents  on  an 
/-plane  (Gill,  1982).  This  radius  is  ~30  km  for  the  model  domain,  as  calculated  by  the 
method  of  Feliks  (1985).  In  addition,  the  equatorward  coastal  jet  and  poleward 
undercurrent  development  in  Experiment  1  agree  well  with  the  /-plane  wind  forcing 
results  of  McCreary  et  al.  (1987)  and  McCreary  (1981). 

Upwelling  along  the  coast  near  the  surface  occurs  during  the  classic  upwelling 
season,  but  offshore  flow  of  the  cooler  upwelled  surface  waters  is  limited,  remaining 
inshore  of  the  intense,  non-meandering,  coastal  equatorward  surface  flow.  The  intensity 
of  the  equatorward  flow  at  depth  results  in  year-round  upwelling  of  deep  waters  offshore, 
on  the  eastern  side  of  the  jet  axis.  A  subsequent  “cold-side”  develops,  encompassing  all 
waters  east  of  the  equatorward  jet. 

As  the  seasons  change,  different  oceanic  responses  do  not  occur  after  the  initial 
stage  of  the  spin-up  phase  due  to  the  continuing  strength  of  the  equatorward  jet.  Surface 
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poleward  flow  in  the  northern  end  of  the  model  domain  does  not  return  in  winter,  as 
expected  in  response  to  both  poleward  wind  stress  in  the  region  and  geostrophic  inflow 
toward  the  coast  from  the  ocean  interior  due  to  the  offshore  pressure  gradient. 

B.  EXPERIMENT  2  -  ROLE  OF  BETA-PLANE 

In  Experiment  2,  the  model  was  run  on  a  (3-plane  with  the  same  realistic  irregular 
coastline,  seasonal  climatological  winds,  thermohaline  forcing  at  the  western  boundary, 
and  annual  salinities  for  MO  as  in  Experiment  1.  The  results  of  Experiment  2  isolate  the 
important  role  that  the  beta  effect  plays,  in  the  presence  of  seasonal  climatological  winds 
and  thermohaline  gradients,  in  setting-up  more  realistic  coastal  currents  that  subsequently 
become  unstable. 

1.  Spin-up  Phase 

Different  oceanic  responses  occur  depending  on  the  season.  In  winter,  the  large 
pressure  gradient  and  poleward  wind  stress  at  the  poleward  end  of  the  model  domain 
result  in  onshore  geostrophic  and  Ekman  flow.  As  the  flow  approaches  the  coast,  it  turns 
and  forms  a  poleward  surface  current.  A  typical  cross-section  of  the  currents  in  the 
poleward  end  of  the  model  domain  (e.g.  Figure  14)  shows  that  the  poleward  surface 
current  exists  within  100  km  of  the  coast  and  has  a  subsurface  core  velocity  of  ~  6  cm/s. 
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When  the  Azores  High  migrates  north  in  the  spring,  the  intensifying  winds  cause 
equatorward  flow  over  the  entire  coastal  domain.  The  weaker  pressure  gradient  and 
stronger  equatorward  wind  stress  at  the  equatorward  end  of  the  model  domain  result  in 
weak  onshore  geostrophic  flow  and  strong  offshore  Ekman  flow. 

As  expected,  upwelling  first  occurs  in  the  south  where  the  stronger  winds  exist 
(not  shown).  In  addition  to  the  west  coast  of  Morocco,  upwelling  is  enhanced  near 
prominent  capes  to  the  north  on  the  west  coast  (e.g.,  off  Cabo  da  Roca  in  Figure  15).  In 
Experiment  1,  the  region  of  enhanced  upwelling  was  confined  to  the  south  off  Morocco 
(e.g.,  Figure  5d).  A  comparison  of  the  temperatures  of  the  coldest  coastal  waters  during 
the  peak  upwelling  season  shows  that  the  upwelled  waters  in  Experiment  2  are  ~1-2°K 
cooler  than  in  the  /-plane  experiment  (Experiment  1)  (e.g.,  compare  Figures  15  and  5d). 
In  addition,  the  offshore  extent  of  the  upwelled  waters  at  the  surface  is  much  wider  in 
Experiment  2. 

A  typical  cross-section  (e.g.,  Figure  16)  of  the  currents  shows  that  there  is  a 
surface  coastal  equatorward  jet  that  extends  to  -180  m  depth  near  the  coast  and  -  800  m 
depth  offshore.  The  surface  current  is  within  -100  km  of  the  coast  and  has  core 
velocities  of  -20-25  cm/s.  Below  the  equatorward  current  is  a  subsurface  poleward 
undercurrent  within  -50  km  of  the  coast  which  has  a  core  velocity  of  -5  cm/s  at  1200m. 
In  contrast  to  the  results  on  an  /-plane  in  Experiment  1,  the  core  of  the  poleward 
undercurrent  is  shallower,  while  the  core  of  the  equatorward  current  has  weaker 
velocities  (e.g.,  compare  Figures  16  and  6a). 
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In  Experiment  1,  there  was  no  development  of  filaments  and,  except  in  the  Gulf  of 
Cadiz,  there  was  also  no  mesoscale  activity.  In  contrast.  Experiment  2  shows  the 
formation  of  meanders,  eddies,  and  filaments  developing  along  the  realistic  irregular 
coastline  (e.g.,  see  the  presence  of  meanders  and  eddies  in  Figure  17).  The  meanders 
preferentially  develop  in  the  vicinity  of  capes  which  subsequently  intensify  to  develop 
predominantly  cold  cyclonic  core  eddies,  which  in  time  coalesce  with  other  cyclonic 
eddies  to  form  larger  cyclonic  eddies  on  the  order  of  ~  100-300  km  diameter. 

In  the  fall,  near  the  end  of  the  upwelling  season,  (e.g.,  Figure  17),  the  equatorward 
surface  flow  has  taken  the  form  of  a  meandering  jet  embedded  with  predominantly 
cyclonic  eddies.  The  meandering  jet  generally  flows  along  the  west  coasts  of  Iberia  and 
Morocco  and  traverses  the  perimeter  of  the  Gulf  of  Cadiz.  While  the  equatorward  jet  in 
Experiment  1  bisected  the  Gulf  of  Cadiz  in  an  anticylonic  manner  to  the  south  and  was 
associated  with  the  development  of  a  cyclonic  eddy  in  the  northern  end  of  the  Gulf  of 
Cadiz,  the  meandering  jet  in  Experiment  2  induces  anticyclonic  flow  at  the  surface  over 
the  entire  Gulf  of  Cadiz  (e.g.,  compare  Figure  17  of  Experiment  2  with  Figure  5b  of 
Experiment  1). 

Of  particular  note  is  the  development  of  an  anticyclonic  eddy  at  ~1200  m  depth 
off  the  southwestern  coast  of  Portugal  near  the  end  of  year  1  (e.g.,  Figure  18a).  A  cross- 
section  of  salinity  at  37.50°N,  in  the  vicinity  of  the  anticyclonic  eddy  (e.g..  Figure  18b) 
shows  a  salty  core  of  ~  35.73  psu  at  -1200  m  depth,  which  is  associated  with  the 
subsurface  poleward  spreading  of  MO.  The  deep  origin,  salty  signature,  depth  and 
anticyclonic  rotation  of  the  eddy  west  of  Cabo  de  Sao  Vincente,  is  consistent  with 
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observations  of  Meddies  in  this  region  (e.g.,  Richardson  and  Tychensky,  1998). 
Maintaining  its  salty  property,  the  anticyclonic  eddy  (Meddy)  subsequently  propagates  to 
the  west  away  from  the  coast  (not  shown)  with  speeds  of  ~2  km/day,  consistent  with 
available  observations  (e.g.,  Richardson  and  Tychensky,  1998).  The  development  and 
subsequent  westward  propagation  of  this  unique  NCCS  feature  was  not  observed  in  the 
results  on  an /-plane  (Experiment  1). 

2.  Quasi-equilibrium  Phase 

Longer  run  times  (~3  years)  of  the  model  simulation  show  that  many  of  the 
features  associated  with  a  classical  eastern  boundary  current  recur  seasonally.  For 
example,  a  surface  poleward  current  develops  during  the  winter  in  the  northern  end  of  the 
model  domain  (e.g.,  Figure  19a),  upwelling  of  cooler  waters  along  the  coast  extends 
offshore  in  the  form  of  cyclonic  eddies  and  narrow  filaments  (e.g.,  Figure  19b),  and  a 
subsurface  poleward  current  beneath  the  meandering  equatorward  jet  is  maintained  (e.g.. 
Figure  19c).  In  addition,  typical  cross-sections  of  the  surface  equatorward  and  subsurface 
poleward  currents  along  the  coast  (e.g..  Figures  19c  and  19d)  show  considerable 
horizontal  and  vertical  shear  in  the  upper  layer  currents. 

To  determine  the  type  of  instability  that  could  generate  the  meander,  eddy,  and 
filament  features,  more  detailed  analysis  was  performed.  Barotropic  instability  can  result 
from  horizontal  shear  in  the  currents,  while  baroclinic  instability  can  result  from  vertical 
shear  in  the  currents.  Horizontal  maps  of  the  upper  layer  mean  kinetic  energy  (MKE)  and 
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eddy  kinetic  energy  (EKE)  are  time  averaged  for  each  month  during  the  upwelling  season 
and  are  used  to  indicate  locations  of  barotropic  and  baroclinic  instability.  Typical  MKE 
and  EKE  results  during  the  upwelling  season  (e.g.,  Figures  20a  and  20b)  show  that  high 
values  of  MKE  and  EKE  are  found  all  along  the  coastal  and  offshore  axes  of  the 
meandering  jet  in  the  regions  of  eddy  and  filament  development,  (e.g.,  compare  Figures 
20a  and  20b  to  Figure  19b).  A  comparison  of  Figures  20a  and  20b  shows  that  maximum 
values  of  MKE  and  EKE  occur  in  the  same  region.  This  is  consistent  with  the  results  of 
Batteen  (1997),  which  showed  that  eddies  are  generated  from  instabilities  of  the  mean 
equatorward  current  and  the  poleward  undercurrent  via  barotropic  and/or  baroclinic 
instability  processes.  Similar  horizontal  maps  of  energetics  for  Experiment  1  (Figures 
20c  and  20d)  show  that  high  values  of  MKE  or  EKE  are  confined  inshore  of  ~12  W  due 
to  the  lack  of  offshore  movement  which  the  beta  effect  produces.  This  suggests  that 
limited  or  no  barotropic  and/or  baroclinic  instabilities  can  exist  when  the  model  is  run  on 
an /-plane. 

During  year  three,  the  Meddies  generated  earlier  in  the  model  run  continue  to 
propagate  to  the  west  and  south  at  ~  2  km/day  (not  shown).  Figure  21a  shows  three 
Meddies  labeled  as  “A”,  “B”  and  “C”.  The  Meddies  have  deep-water  origins  at  -1200  m 
depth  off  southwest  Portugal.  Cross-sections  of  meridional  velocity  (Figure  21b)  and 
temperature  (Figure  21c)  associated  with  Meddy  “A”  illustrate  the  anticyclonic  rotation 
of  the  Meddy.  The  corresponding  cross-section  of  salinity  (Figure  21d)  shows  the 
relatively  salty  signature  of  the  Meddy  at  -1200  m  depth  which  it  retains  as  it  propagates 
westward  (e.g.,  see  the  high  salinity  cores  at  -11.5°W  and  between  15°W  and  16°W  in 
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Figure  2 Id).  Unlike  Experiment  2,  no  Meddy  development  off  Cabo  de  Sao  Vincente 
was  observed  in  Experiment  1.  Also,  no  westward  propagation  of  any  anticyclonic  or 
cyclonic  eddies  occurred  in  Experiment  1. 

3.  Discussion 

To  help  illustrate  the  importance  of  Coriolis  parameterization  (specifically  the 
role  of  planetary  beta),  in  the  presence  of  alongshore  and  cross-shore  varying  wind  fields 
and  thermohaline  gradients,  in  developing  realistic  EBC  and  unique  NCCS  features,  the 
following  summary  of  Experiment  2  results  is  provided: 

-  Responsiveness  to  seasonal  oceanic  changes  (e.g.,  development  of  a  poleward 

surface  current). 

-  Realistic  locations,  depths,  and  velocities  of  coastal  currents. 

-  Realistic  offshore  extent  of  upwelling  patterns  with  enhanced  upwelling  in  the 

vicinity  of  capes. 

-  Upper  layer  barotropic  and  baroclinic  instability. 

-  Typical  EBC  mesoscale  activity  (i.e.  eddies,  meandering  jet,  filaments). 

-  Poleward  movement  of  MO,  Meddy  generation,  and  subsurface  westward 

propagation  of  Meddies. 

Since  none  of  these  features  were  present  in  the  results  of  Experiment  1,  it  can  be 
concluded  that  their  existence  is  related  to  the  beta  effect. 
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The  beta  effect  allows  the  existence  of  freely  propagating  planetary  waves,  i.e., 
Rossby  waves  (Gill,  1982).  The  offshore  propagation  of  these  waves  contributes  to  the 
generation  of  an  alongshore  pressure  gradient  field,  which  aids  the  development  of  a 
realistic  undercurrent  along  the  eastern  boundary.  In  addition,  the  offshore  extent  of  the 
surface  equatorward  current  (e.g.,  Figure  19c)  was  due  to  the  beta  effect.  This  is 
consistent  with  the  results  of  McCreary  et  al.,  1987,  which  showed  that,  due  to  the  beta 
effect,  the  surface  coastal  jet  does  not  necessarily  have  to  be  confined  to  within  a  Rossby 
radius  of  deformation  of  the  coast.  As  a  result,  the  beta  effect  changes  both  the  vertical 
and  horizontal  structures  of  the  currents. 

In  particular,  as  the  core  of  the  undercurrent  (which  is  initially  beneath  the 
equatorward  jet)  intensifies,  it  shoals  and  displaces  the  equatorward  jet  offshore  (e.g., 
compare  Figures  16  and  19c).  Strong  vertical  and  horizontal  shears  in  the  upper  layers 
between  the  equatorward  jet  and  the  undercurrent  result  (e.g.  see  Figures  19c  and  19d), 
and  meanders,  upwelling  filaments,  and  eddies  subsequently  develop  in  the  baroclinically 
and  barotropically  unstable  coastal  environment.  Note  that  neither  vertical  nor  horizontal 
shears  existed  between  the  deep  and  intense  equatorward  jet  and  the  even  deeper 
poleward  undercurrent  in  Experiment  1  (/-plane).  As  a  result,  the  equatorward  jet  was 
unmodified  by  the  poleward  subsurface  flow,  and  no  meanders,  filaments  nor  eddies 
developed  in  Experiment  1 . 

The  beta  effect  also  played  an  important  role  in  the  unique  development  and 
westward  propagation  of  Meddies  from  the  MO  in  the  NCCS.  It  is  suggested  that  the 
Meddies  are  generated  by  the  basic  instability  of  the  equatorward  coastal  jet  and  the 
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poleward  undercurrent.  The  major  difference  between  Meddies  and  the  previously 
discussed  eddies  is  that  the  Meddies  are  generated  at  depth  rather  than  near  the  surface. 
The  relatively  salty  MO  core  of  the  undercurrent  below  -500  m  depth,  which  makes 
water  in  the  core  saltier  than  water  above  or  below  it,  could  lead  to  a  local  baroclinic 
instability  process,  resulting  in  the  generation  of  Meddies.  This  scenario  is  consistent 
with  the  results  of  Kase  et  al.,  (1989),  who  observed  eddies  in  the  MO  off  Iberia  similar 
to  the  Meddy  off  Cabo  de  Sao  Vincente  in  Figure  18a  and  postulated  that  baroclinic 
instabilities  generated  by  the  undercurrent  could  lead  to  the  generation  of  Meddies. 

The  development  of  the  Meddy  off  Cabo  de  Sao  Vincente  in  Figure  18a  is  also 
consistent  with  another  Meddy  generation  theory.  Pichevin  and  Nof  (1996)  suggested 
that  a  Meddy  generation  mechanism  can  occur  when  the  MO  curves  back  on  itself  west 
of  Cabo  de  Sao  Vincente  producing  a  flow-force  that  cannot  be  balanced  without  growth 
and  shedding  of  eddies.  The  beta  effect  plays  an  important  role  in  the  Meddy  generation 
model  of  Pichevin  and  Nof  (1996)  because  it  arrests  the  growth  of  Meddies  and  removes 
them  from  their  generation  region.  (In  contrast, /-plane  experiments  by  Pichevin  and  Nof 
(1996)  develop  a  constantly  growing  coastal  eddy  at  the  tip  of  Cabo  de  Sao  Vincente.) 

In  Experiment  2,  the  beta  effect  acted  as  a  necessary  force  in  promoting  the 
generation  and  growth  of  Meddies  and  causing  their  removal  from  the  source  region. 
Generation  and  growth  were  promoted  by  local  baroclinic  instabilities  in  the  MO  and  by 
the  MO  curving  back  on  itself  in  a  region  of  relatively  weak  coastal  currents.  Removal 
from  the  Cabo  de  Sao  Vincente  generation  region  was  a  result  of  the  free  offshore 
propagation  of  planetary  waves,  in  the  vicinity  of  the  relatively  benign  coastal  currents. 
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In  contrast  to  both  Experiment  2  and  the /-plane  experiment  of  Pichevin  and  Nof  (1996), 
Meddy  generation  was  not  simulated  in  Experiment  1  due  to  the  domination  of  the 
intense  equatorward  jet  through  the  potential  Meddy  generation  region.  In  addition,  the 
waters  with  a  salty  MO  signature  in  Experiment  1,  instead  of  freely  propagating  offshore, 
were  confined  to  the  Gulf  of  Cadiz  and  the  narrow  region  between  the  equatorward  jet 
and  west  facing  coastlines  (e.g.,  see  Figure  12). 

To  more  realistically  investigate  the  effects  of  MO  on  the  NCCS,  two  additional 
experiments  (Experiments  3  and  4)  were  run  using  seasonal  values  of  salinity  for  the  MO 
(see  Figure  4).  In  addition,  the  Strait  of  Gibraltar  was  opened  (see  Figure  lb)  to  allow  for 
the  free  exchange  of  water  with  the  Mediterranean  Sea. 

Results  for  the  /-plane  case  (Experiment  3),  showed  similar  effects  as  the  /-plane 
case  (Experiment  1),  i.e.,  the  intense  equatorward  jet  dominated  the  entire  NCCS  region 
and  stifled  the  development  of  upper  layer  mesoscale  activity  and  Meddies.  Experiment 
4,  run  on  a  (5-plane,  also  produced  results  similar  to  the  3-plane  case  (Experiment  2)  in 
the  upper  layers.  Since  it  more  realistically  incorporated  the  influences  of  MO  at  depth 
along  with  the  realistic  opening  of  the  Strait  of  Gibraltar,  it  will  be  compared  with 


observations. 


C.  COMPARISON  WITH  OBSERVATIONS 


It  is  useful  to  qualitatively  compare  Experiment  4  model  simulation  results  with 
observational  data  because,  of  all  the  experiments,  it  most  realistically  depicts  the 
classical  and  unique  features  of  the  NCCS.  Since  this  experiment  is  an  idealized  process- 
oriented  study  and  not  a  model  hindcast,  we  cannot  make  direct  comparisons  with  data; 
however  we  can  investigate  whether  the  phenomenological  model  behavior  is 
qualitatively  similar  to  observational  data  in  the  NCCS.  Overall,  the  model  behavior  of 
this  experiment’s  complex  flow  regime  highlights  the  mesoscale  features  of  the  NCCS 
with  relatively  close  similarities  to  field  observations  in  both  the  upper  and  lower  layers. 

In  previous  studies,  major  EBC  features  produced  by  the  PE  model  were  similar 
to  observations  of  mesoscale  activity  in  the  upper  layers  (e.g.  see  Batteen,  2000).  In 
similar  fashion,  Experiment  4  captures  these  classical  EBC  features.  The  mean 
equatorward  surface  flow  of -10-30  cm/s  throughout  the  model  domain  is  consistent  with 
Wooster  and  Reid  (1963)  and  observations  by  DMAH/TC  (1988).  The  poleward  flow 
along  Iberia  developed  in  the  model  is  in  conformity  with  hydrographic  survey  data  from 
1986  (Haynes  and  Barton,  1990),  and  the  poleward  surface  flow  depicted  in  the  model  is 
in  accord  with  infrared  satellite  images  and  in  situ  hydrographic  data  from  Frouin  et  al., 
(1990). 

The  upwelling  patterns  and  cold  filaments  in  Experiment  4  were  consistent  with 
high  resolution  satellite  SST  observations  of  these  features  off  Cape  Ghir,  Cape 
Beddouzza  and  southwestern  Portugal  (e.g.,  Van  Camp  et  al,  1991),  and  off  western 
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Portugal  (Foulkard  et  al,  1997).  In  addition,  the  results  of  Experiment  4  show 
termination  of  filaments  in  cyclonic  and  anticyclonic  pairs  of  eddies  as  observed  in 
satellite  observations  (e.g.,  Van  Camp  et  al.,  1991  and  Hagen  et  al.,  1996). 

In  the  lower  layers,  the  behavior  of  MO  as  it  flows  along  the  southern  coast  of 
Iberia  in  the  Gulf  of  Cadiz  and  around  Cabo  de  Sao  Vincente  in  two  distinct  cores  is  well 
documented  (e.g.  Ambar  and  Howe,  1979;  Ambar,  1982;  Fiuza,  1982;  Zenk  and  Armi, 
1990;  Iorga  and  Lozier,  1999).  CTD  stations  near  Cabo  de  Sao  Vincente  by  Daniault  et 
al.,  (1994)  show  upper  core  and  lower  core  high  salinity  maxima  and  support  this  two 
core  phenomenon.  Typical  cross-sections  of  salinity  in  Experiment  4  near  Cabo  de  Sao 
Vincente  show  similar  high  salinity  values  associated  with  the  poleward  and  westward 
spreading  of  MO  at  -800  and  -1200  m  depth  (e.g.,  Figure  22).  In  addition,  typical  cross- 
sections  of  salinity  throughout  the  model  domain  are  relatively  consistent  with  the 
assembled  database  of  the  North  Atlantic  Ocean  from  the  National  Ocean  Data  Center  by 
Lozier  et  al.  (1995)  that  traces  the  salinity  signature  of  MO  throughout  the  NCCS  region. 

In  Experiment  4,  the  salty  characteristics  and  structure  of  Meddy  “A”  in  Figures 
23a  through  23e  are  consistent  with  recent  observations  of  Meddies  in  the  Canary  Basin 
(e.g.  Prater  and  Sanford,  1994;  Richardson  and  Tychensky,  1998;  Stammer  et  al.,  1991). 
Rotation  speeds  of  between  -10-20  cm/s  (e.g..  Figures  23a,  23b,  and  23e)  and  the 
elliptical  orientation  along  a  northeast  to  southwest  line  in  the  lower  part  of  the  Meddy 
(e.g.,  Figure  23e)  are  consistent  with  observations  by  Prater  and  Sanford  (1994). 
Propagation  speeds  and  propagation  paths  of  Meddies  in  Experiment  4  agree  with 
observations  by  Richardson  and  Tychensky  (1998).  Anticyclonic  motion  near  the  surface 
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above  Meddy  “A”  in  Experiment  4  (not  shown)  is  in  agreement  with  observations  of 
Meddies  made  by  Stammer  et  al  (1991)  which  compared  a  Geosat  signal  with  the 
hydrographic  survey  data  of  Kase  et  al.  (1989)  in  the  NCCS  region. 

In  the  model  simulations,  the  influences  of  the  beta  effect  and  the  associated 
offshore  movement  of  planetary  Rossby  waves  resulted  in  accurate  depiction  of  observed 
mesoscale  features  of  the  NCCS.  Using  Geosat  altimeter  data,  Tokmakian  and 
Chancellor  (1993)  observed  baroclinic  Rossby  waves  between  ~30°N  and  40°N  in  the 
vicinity  of  the  NCCS  that  appeared  to  originate  off  the  African  coast. 
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IV.  SUMMARY 


The  objective  of  this  study  was  to  use  a  multilevel  PE  model  to  investigate  the 
role  of  planetary  beta  on  classical  as  well  as  unique  features  in  the  predominantly  wind- 
driven  NCCS.  Toward  this  end,  four  numerical  experiments  were  run  with  varying 
Coriolis  parameterizations  (/"-plane  or  3-plane).  The  first  two  experiments  used  a  closed 
boundary  and  annual  salinity  forcing  for  MO.  The  latter  two  experiments  used  an  open 
Mediterranean  Sea  at  the  Strait  of  Gibraltar  and  seasonal  salinity  forcing  for  MO  to 
permit  a  more  accurate  investigation  of  the  role  of  the  beta  effect  on  subsurface  spreading 
of  MO  and  Meddies. 

In  a  similar  fashion  to  classical  EBC  systems,  both  experiments  run  on  a  3-plane 
(Experiments  2  and  4)  produced  an  offshore  surface  equatorward  meandering  jet,  realistic 
surface  and  subsurface  poleward  currents,  upwelling,  meanders,  eddies,  and  filaments.  In 
addition,  these  experiments  depicted  unique  NCCS  features  including  geographical 
separation  of  the  of  the  Gulf  of  Cadiz  region  from  the  west  coast  upwelling  regimes, 
anticyclonic  mesoscale  activity  in  the  upper  layers  of  the  Gulf  of  Cadiz,  poleward 
spreading  of  MO,  and  the  development  and  propagation  of  Meddies  from  the  Cabo  de 
Sao  Vincente  generation  region. 

To  isolate  the  role  of  the  beta  effect,  the  results  of  the /-plane  runs  (Experiments  1 
and  3)  were  compared  to  the  3 -plane  runs  (Experiments  2  and  4).  Results  on  an /-plane 
showed  the  occurrence  of  coastal  upwelling  and  the  eventual  appearance  of  anticyclonic 
mesoscale  activity  in  the  Gulf  of  Cadiz,  but  the  continuously  strengthening  equatorward 
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jet  prevented  other  classical  mesoscale  features  associated  with  upwelling  systems  from 
occurring.  Specifically,  upwelling  was  confined  to  the  immediate  coastal  regions,  a 
seasonal  surface  poleward  current  did  not  occur  in  the  north,  the  subsurface  poleward 
flow  was  forced  to  great  depths,  barotropic  and  baroclinic  instabilities  were  relatively 
non-existent,  and  meanders,  filaments,  and  eddies  did  not  develop.  At  depth,  the  deep 
equatorward  jet  prevented  unique  NCCS  features  from  occurring,  including  the  typical 
westward  and  poleward  spreading  of  MO  and  the  development  of  Meddies.  In  addition, 
cold  upwelled  water  and  salty  MO  were  restricted  to  the  region  east  of  the  equatorward 
jet  and  in  the  Gulf  of  Cadiz. 

The  additional  effects  of  seasonally  forcing  MO  and  opening  the  eastern  boundary 
at  the  Strait  of  Gibraltar  in  Experiment  4  more  realistically  depicted  the  poleward 
spreading  of  MO  and  produced  sub-mesoscale  coherent  vortices  (SCVs)  off  Cabo  de  Sao 
Vincente  with  more  typical  Meddy  characteristics.  Specifically,  Meddies  in  Experiment 
4  were  shed  regularly  from  the  classic  generation  region  (in  numbers  consistent  with 
available  observations),  and  with  rotation  velocities,  high  salinities,  and  temperature 
patterns  consistent  with  various  Meddy  observations. 

Meddies  produced  near  Cabo  de  Sao  Vincente  in  the  (3-plane  experiments 
supported  recent  theories  on  Meddy  generation.  In  particular,  the  generation  mechanisms 
of  Meddies  appeared  to  be  due  to  baroclinic  instabilities  in  the  undercurrent  (saltier  MO 
cores  than  the  water  above  or  below)  as  suggested  by  Kase  et  al.,  1989.  In  addition,  the 
Meddies  also  appeared  to  be  generated  due  to  the  curvature  of  the  MO  back  on  itself  (i.e., 
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the  flow  force  could  not  be  balanced  without  generating  and  shedding  anticyclonic 
eddies),  as  suggested  by  Pichevin  and  Nof  (1996). 

A  comparison  of  currents,  upwelling,  meanders,  eddies,  jets,  filaments,  MO,  and 
Meddies  generated  by  the  model  in  Experiment  4  with  available  observations  showed  that 
the  model  successfully  simulated  the  location,  size,  and  velocity  of  these  features  when 
run  on  a  (3 -plane.  In  particular,  the  3 -plane  results  of  the  complex  upper  layer  and 
subsurface  flow  regimes  of  Experiment  4  showed  that  the  incorporation  of  the  beta  effect 
resulted  in  an  accurate  depiction  of  observed  mesoscale  activity,  and  highlighted  major 
classical  EBC  and  unique  NCOS  mesoscale  features. 
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Figure  1.  The  model  domain  for  the  northern  Canary  Current  System  (NCCS)  is  bounded 
by  30°N  to  42.5°N,  5°W  to  17.5°W.  The  cross-shore  (alongshore)  resolution  is  9  km  (11 
km).  Geographical  locations  and  prominent  features  are  labeled.  Model  domain  in  (a) 
has  a  closed  boundary  along  the  entire  coast,  and  in  (b)  has  been  modified  at  5  grid  points 
at  the  eastern  boundary  to  open  the  Mediterranean  Sea  at  the  Strait  of  Gibraltar. 
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Figure  2.  Climatological  (1980-1989)  ECMWF  winds  in  m/s  for:  (a)  December,  (b) 
February,  (c)  April,  (d)  August,  and  (e)  October.  Maximum  wind  vector  is  8  m/s. 
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Figure  3.  Plot  of  Levitus  climatological  fields  used  as  seasonal  forcing  along  the  western 
boundary  (17.5°W)  at  42.5°N  (top)  and  30°N  (bottom)  for  (a)  temperature  and  (b)  salinity. 


44 


Salinity  (psu)  Salinity  (psu) 


Salinity  Forcing  Along  Western  Boundary  (17.5W) 


Months  (January  -  December) 
(b) 


5 


Salinity  (PSU) 


Days 


Figure  4.  Salinities  for  Mediterranean  Outflow  (MO)  used  as  seasonal  forcing  at  the 
Strait  of  Gibraltar  (from  Levitus  et  al.,  1994). 
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Figure  5a.  Temperature  contours  and  velocity  vectors  at  10  m  depth  for  Experiment  1  at 
day  27.  To  avoid  clutter,  velocity  vectors  are  plotted  every  third  gridpoint  in  the  cross¬ 
shore  and  alongshore  directions.  Contour  interval  is  1°  C.  Maximum  velocity  vector  is 
100  cm/s. 
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Figure  5b.  Velocity  vectors  at  10  m  depth  for  Experiment  1  at  day  267.  To  avoid  clutter, 
velocity  vectors  are  plotted  every  third  gridpoint  in  the  cross-shore  and  alongshore 
directions.  Maximum  velocity  vector  is  100  cm/s. 
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Figure  5c.  Cross-shore  section  of  meridional  velocity  (v)  at  35.3°N  for  Experiment  1  on 
day  330.  Contour  interval  is  10  cm/s.  Equatorward  (poleward)  flow  is  denoted  by 
dashed  (solid)  lines. 
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Figure  5d.  Temperature  contours  and  velocity  vectors  at  10  m  depth  for  Experiment  1  at 
day  216.  To  avoid  clutter,  velocity  vectors  are  plotted  every  third  (fourth)  gridpoint  in 
the  cross-shore  (alongshore)  direction.  Contour  interval  is  1°  C.  Maximum  velocity 
vector  is  100  cm/s. 
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Figure  6a.  Cross-shore  section  of  meridional  velocity  (v)  at  37.3°N  for  Experiment  1  on 
day  231.  Contour  interval  is  10  cm/s  for  equatorward  flow  (dashed  lines)  and  2  cm/s  for 
poleward  flow  (solid  lines). 
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Figure  6b.  Cross-shore  section  of  meridional  velocity  (v)  at  41.5°N  for  Experiment  1  on 
day  357.  Contour  interval  is  10  cm/s  for  equatorward  flow  (dashed  lines,  with  an 
additional  5  cm/s  contour  added)  and  5  cm/s  for  poleward  flow  (solid  lines). 
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Figure  7.  Cross-shore  section  of  salinity  for  Experiment  1  at  37.1°N  on  day  273.  The 
contour  interval  for  salinity  is  0.01  (0.02)  above  (below)  35.69  values. 
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Figure  8a.  Velocity  vectors  at  10  m  depth  for  Experiment  1  on  day  519.  To  avoid 
clutter,  velocity  vectors  are  plotted  every  third  gridpoint  in  the  cross-shore  and 
alongshore  directions.  Maximum  velocity  vector  is  100  cm/s. 
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Figure  8b.  Velocity  vectors  at  10  m  depth  for  Experiment  1  at  day  561.  To  avoid  clutter, 
velocity  vectors  are  plotted  every  third  gridpoint  in  the  cross-shore  and  alongshore 
directions.  Maximum  velocity  vector  is  100  cm/s. 
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Figure  8c.  Velocity  vectors  at  10  m  depth  for  Experiment  1  at  day  942.  To  avoid  clutter, 
velocity  vectors  are  plotted  every  third  gridpoint  in  the  cross-shore  and  alongshore 
directions.  Maximum  velocity  vector  is  100  cm/s. 
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Figure  9.  Temperature  contours  and  velocity  vectors  at  10  m  depth  for  Experiment  1  at 
day  999.  To  avoid  clutter,  velocity  vectors  are  plotted  every  fifth  gridpoint  in  the  cross¬ 
shore  and  alongshore  direction.  Contour  interval  is  1°  C.  Maximum  velocity  vector  is 
100  cm/s. 
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Figure  10a.  Cross-shore  section  of  meridional  velocity  (v)  at  39.3°N  for  Experiment  1  on 
day  762.  Contour  interval  is  10  cm/s  for  equatorward  flow  (dashed  lines)  and  5  cm/s  for 
poleward  flow  (solid  lines). 
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Figure  11a.  Cross-shore  section  of  temperature  for  Experiment  1  at  39.3°N  on  day  828 
The  contour  interval  for  temperature  is  0.5°C. 
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Figure  lib.  Time-averaged  plot  of  temperature  contours  and  velocity  vectors  at  10  m 
depth  for  Experiment  1  in  the  third  year  for  October.  To  avoid  clutter,  velocity  vectors 
are  plotted  every  fifth  gridpoint  in  the  cross-shore  and  alongshore  directions.  Contour 
interval  is  2°  C.  Maximum  velocity  vector  is  100  cm/s. 
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Figure  12.  Time-averaged  plot  of  salinity  contours  and  velocity  vectors  at  600  m  depth 
for  Experiment  1  in  the  third  year  for  October.  To  avoid  clutter,  velocity  vectors  are 
plotted  every  sixth  gridpoint  in  the  cross-shore  and  alongshore  directions.  Contour 
interval  is  0.04.  Maximum  velocity  vector  is  100  cm/s. 
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Figure  13.  Cross-shore  section  of  salinity  for  Experiment  1  at  37.1°N  on  day  1002.  The 
contour  interval  for  salinity  is  0.01  (0.02)  above  (below)  35.69  values. 
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Figure  14.  Cross-shore  section  of  meridional  velocity  (v)  at  40.3°N  for  Experiment  2  on 
day  363.  Contour  interval  is  10  cm/s  for  equatorward  flow  (dashed  lines)  and  2  cm/s  for 
poleward  flow  (solid  lines). 
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Figure  15.  Temperature  contours  and  velocity  vectors  at  10  m  depth  for  Experiment  2  at 
day  216.  To  avoid  clutter,  velocity  vectors  are  plotted  every  third  (fourth)  gridpoint  in 
the  cross-shore  (alongshore)  direction.  Contour  interval  is  1°  C.  Maximum  velocity 
vector  is  100  cm/s. 
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Figure  16.  Cross-shore  section  of  meridional  velocity  (v)  at  37.3  °N  for  Experiment  2  on 
day  210.  Contour  interval  is  5  cm/s  for  equatorward  flow  (dashed  lines)  and  2.5  cm/s  for 
poleward  flow  (solid  lines). 
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Figure  18b.  Cross-shore  section  of  salinity  for  Experiment  2  at  37.5°N  on  day  273.  The 
contour  interval  for  salinity  is  0.02  (0.03)  above  (below)  35.63  values. 
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Figure  19a.  Cross-shore  section  of  meridional  velocity  (v)  at  41.5°N  for  Experiment  2  on 
day  1089.  Contour  interval  is  5  cm/s  for  equatorward  flow  (dashed  lines)  and  for 
poleward  flow  (solid  lines). 
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Figure  19c.  Cross-shore  section  of  meridional  velocity  (v)  at  39.3°N  for  Experiment  2  on 
day  909.  Contour  interval  is  10  cm/s  for  equatorward  flow  (dashed  lines)  and  2.5  cm/s 
for  poleward  flow  (solid  lines). 
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Figure  20a.  Horizontal  map  at  10  m  depth  of  average  mean  kinetic  energy  (MKE)  for 
Experiment  2  in  the  third  year  of  the  model  simulation,  time  averaged  for  July  through 
September.  Contour  interval  is  200  (cm/s)2 . 
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Figure  20b.  Horizontal  map  at  10  m  depth  of  average  eddy  kinetic  energy  (EKE)  for 
Experiment  2  in  the  third  year  of  the  model  simulation,  time  averaged  for  July  through 
September.  Contour  interval  is  200  (cm/s)2 . 
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Figure  20c.  Horizontal  map  at  10  m  depth  of  average  mean  kinetic  energy  (MKE)  for 
Experiment  1  in  the  third  year  of  the  model  simulation,  time  averaged  for  July  through 
September.  Contour  interval  is  1000  (cm/s)2. 
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Figure  20d.  Horizontal  map  at  10  m  depth  of  average  eddy  kinetic  energy  (EKE)  for 
Experiment  1  in  the  third  year  of  the  model  simulation,  time  averaged  for  July  through 
September.  Contour  interval  is  200  (cm/s)2 . 
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Figure  21a. .  Velocity  vectors  at  1226  m  depth  for  Experiment  2  at  day  1059.  Velocity 
vectors  are  plotted  every  second  gridpoint  in  the  cross-shore  and  alongshore  directions. 
Maximum  velocity  vector  is  20  cm/s.  Meddies  are  labeled  “A”,  “B”,  and  “C”. 
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Figure  21b.  Cross-shore  section  of  meridional  velocity  (v)  at  35.5°N  for  Experiment  2 
on  day  1059.  Contour  interval  is  5  cm/s  for  equatorward  flow  (dashed  lines)  and 
poleward  flow  (solid  lines). 
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Figure  22.  Cross-shore  section  of  salinity  for  Experiment  4  at  37.1°N  on  day  273.  The 
contour  interval  for  salinity  is  0.01  (0.02)  above  (below)  35.63  (35.60)  values. 
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Figure  23a. .  Velocity  vectors  at  1226  m  depth  for  Experiment  4  at  day  1068.  Velocity 
vectors  are  plotted  every  second  gridpoint  in  the  cross-shore  and  alongshore  directions. 
Maximum  velocity  vector  is  20  cm/s.  Meddy  “A”  is  labeled. 
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Figure  23c.  Cross-shore  section  of  temperature  for  Experiment  4  at  37.7°N  on  day  1068. 
The  contour  interval  for  temperature  is  0.2  °C. 
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Figure  23e.  Temperature  contours  and  velocity  vectors  for  Experiment  4  at  1226  m  depth 
on  day  1068.  The  contour  interval  is  0.02°C.  The  maximum  current  velocity  is  20  cm/s. 
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Table  1.  Summary  of  Batteen  (1989)  specific  experimental  design. 


Cases 

Coastline 

Coriolis 

Parameterization 

Type  of  Wind 
Stress  (steady  wind) 

1 

/-plane 

constant 

2 

■EB39H 

3 -plane 

constant 

3 

/-plane 

4 

■KBS9HI 

3-plane 

Table  2.  Summary  of  Batteen  (1997)  specific  experimental  design. 


Cases 

Coastline 

Coriolis 

Parameterization 

Type  of  Wind  Stress 
(steady  wind) 

1 

3 -plane 

2a 

2b 

/-plane 

2c 

■EBD9H 

■bsh 

constant 

2d 

/-plane 

constant  (interior  only) 

3 

M 

3-plane 

Table  3.  Summary  of  Batteen(2000)  specific  experimental  design. 


Cases 

Coastline 

Coriolis 

Parameterization 

Annual  Mediterranean  Outflow 
and  Thermohaline  Forcing 

1 

3-plane 

NO  | 

2 

3 -plane 

NO  i 

3 

3 -plane 

NO 

4 

3 -plane 

YES 

Table  4.  Summary  of  specific  experimental  design  in  this  study. 


Cases 

Coastline 

Coriolis 

Parameterization 

Thermohaline 

Gradients 

Annual  or 
Seasonal  MO 

1 

I^^KSBaNHI 

YES 

Annual 

2 

usaimnsHsira 

3 -plane 

YES 

Annual 

3 

mmsmmm 

YES 

Seasonal 

4 

3 -plane 

YES 

Seasonal 
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Table  5.  Values  of  constants  used  in  the  model 


I  Constant 

Value 

Definition 

To 

278.2°K 

Constant  Reference  Temperature 

So 

34.7 

Constant  Reference  Salinity 

Po 

Density  of  Sea  Water  At  T0  and  S0 

a 

2.4  x  IQ-4  (0K)4 

Thermal  Expansion  Coefficient 

P 

Saline  Expansion  Coefficient 

K 

10 

Number  of  Levels  In  Vertical 

Ax 

Cross-Shore  Grid  Spacing 

Ay 

mmBmm 

Alongshore  Grid  Spacing 

H 

Total  Ocean  Depth 

At 

800  s 

Time  Step 

fo 

Mean  Coriolis  Parameter 

8 

BM 

Acceleration  of  Gravity 

A-m 

EGQB3SH 

Biharmonic  Momentum  Diffusion  Coefficient 

Ah 

1 

Biharmonic  Heat  Diffusion  Coefficient 

0.5  cm2  s‘A 

Vertical  Eddy  Viscosity 

kh 

0.5  cm2  s'1 

Vertical  Eddy  Conductivity 

90 


Table  6:  Annual  salinity  profiles  for  MO 


Layer 

(Depth  in  meters) 

Salinity  (psu) 

1 

10 

36.07 

2 

30 

36.05 

3 

75 

36.05 

150 

36.05 

5 

250 

36.05 

6 

400 

35.75 

7 

600 

35.75 

8 

1226 

36.25 
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